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The Strange Quark Mass, 0^5, and the Chiral Limit Electroweak Penguin 
K ^ TTTT Matrix Elements From Hadronic r Decay Data 
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Hadronic r decay data provides access to the light quark vector (V) and axial vector (A) spectral hinctions. 
This makes possible investigations of the dynamics of QCD at intermediate scales and improved determinations 
of certain QCD/Standard Model parameters. We discuss three such applications: (1) the investigation of the 
nature of duality violation in QCD at intermediate scales (and its relation to the determination of as); (2) the 
extraction of from flavor-breaking differences of ud and us data; and (3) the determination of the dimension 
D — 6 term in the OPE for the flavor ud V-A correlator difference. The latter is relevant to the evaluation of 
the chiral limit values of the K tttt electroweak penguin matrix elements, and hence to our understanding of 
expectations for e'/e in the Standard Model. 
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1. INTRODUCTION/BACKGROUND 



As is well known, the ratio, R 



V/A:ij 



of the fla- 



vor ij = ud, us vector (V) or axial vector (A) 
current-mediated hadronic t decay rate to the 
corresponding electron decay rate. 



R 



V/A-tj 



T[t Ur hadronsv/A;ij (7)] 



r[r- 



(1) 



(with (7) denoting extra photons and/or lepton 
pairs) can be written as a sum of weighted inte- 
grals over the spin J = and 1 hadronic spectral 
functions . Explicitly, 



R 



[127r2|U,,|2W] JO 



dVr (1 - Vr) 



\^ + ^yr)p^v;Ah{s)~2yJ^U,^{s) 
dRv/A-ij (s) 



ds ■ 



ds 



(2) 



where y-r = s/rn^, Vij is the flavor ij CKM ma- 
trix element. Sew is an electroweak correction, 
and Py'/A-ijis), with (J) the spin of the hadronic 
system, is the spectral function of the correspond- 
ing spin (J) part of the flavor ij V/A correlator, 
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^v/A■,^,i^)^ The n|,°;2^^.(s) are defined by 
I j d^xe^^-^Q\T(j'^^^{x)ry/^{Q)) |0) 



V/A 



where Jy^^ are the standard V and A currents. 
For a given channel, the kinematically-weighted 
linear combination of (J) = (0 -I- 1) and (0) spec- 
tral functions appearing in Eq. (H) can thus be 
extracted from the corresponding bin-by-bin ex- 
perimental decay distribution. 

For the V channel, Py^^j is proportional to 
(rrii — TOj)^, and hence numerically negligible for 
ij = ud. For the A channel, p^2\j is saturated by 
the flavor ij Goldstone boson pole in the chiral 
limit and has non-Goldstone-boson contributions 
proportional to (m^ -I- rrijY. Apart from the tt 

pole contribution, p^^.^j^ is, thus, also numerically 

negligible, mg is not sufficiently small that Pv\s 

and the non-_fi'-pole contributions to p''^^^ can 
be safely neglected. The spin separation of the 

2p(0)(s) + = p(0+l)(s). This combination, and 

,r correlators, nj*?^/' . = 

having no kinematic sin- 



also sp''')(s), correspond to scalar correlators, bl^^"!^^'^. 
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us experimental data is straightforward for the 
K and K* contributions, but not yet available 
for the data above the K*. Flavor ud V/A sep- 
aration for states consisting only of pions can be 
accomplished using G-parity. 

Cauchy's theorem, together with analytic- 
ity, implies that correlators 11 having no kine- 
matic singularies satisfy finite energy sum rules 
(FESR's). For weight functions, w{s), analytic in 
the region \s\ < S oi the complex-s plane, and 
any sq < S, these have the form 

so _i r 

dsp{s)w{s) — - — : f dsll{s) w{s) .(4) 
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For So large enough that 11 on the RHS of Eq. 
may be approximated by its OPE, one obtains a 
relation between spectral data and OPE param- 
eters. The well-known predictions for the inclu- 
sive flavor ij = ud, us hadronic r decay widths 
in terms of the parameters occurring in the OPE 
of the corresponding V and A correlators (domi- 
nantly a^) 



dy (1 - y) 



(i + 2y)<+^^(s)-2ynV7^^^^.(s)j (5) 

are examples of such a relation. The modified ver- 
sion of this relation obtained by multiplying the 
integrands in both the spectral and OPE repre- 
sentations of Rv/A-ij by the factor (1 — yr)''y^ is 
called "the (fc, m) spectral weight sum rule" |]. 
Such sum rules have the advantage that a spin 
separation of the ij = us spectral data is not re- 
quired in order to evaluate their spectral sides. 

The spectral integrals required for various V 
and A FESR's can be evaluated, using hadronic 
r decay data, for any sq < ra^. So long as there 
exists a window of sq values below for which 
the breakdown of the OPE representation of the 
relevant correlator ("duality violation") is negli- 
gible, the values of basic QCD parameters, such 
as as and , appearing on the OPE side of the 
FESR relation can be determined using spectral 
integral data. The existence of such a window of 
So values is, of course, crucial to the reliability 
of such a determination. We will concentrate on 
three analyses of this type below. 



(0) 



The rest of the paper is organized as follows. 
In Section 2 we discuss the issue of duality vio- 
lation, and the closely related issue of the deter- 
mination of as, using the very precise ud V-l-A 
spectral data. In Section 3, we discuss some com- 
plications, and the current status, of attempts to 
extract ttt-s from the flavor-breaking difference of 
ud and us V+A correlators. In Section 4, we 
discuss some recent work on the determination 
of the dimension 6 term in the OPE of the ud 
V-A correlator. This quantity bears a special re- 
lation to the chiral limit values of the K tttt 
electroweak penguin (EWP) operator matrix ele- 
ments, and hence to expectations for the value of 
e'/e in the Standard Model. Finally, in Section 5, 
we consider what improvements in these determi- 
nations may be possible in the near future. 

2. DUALITY VIOLATION AND THE 
EXTRACTION OF 

For very large sq, the OPE is expected to pro- 
vide a reliable representation of hadronic correla- 
tors over the entire circle |s| = sq in the complex 
s-plane. For such sq, the corresponding spectral 
function, p(so), will also be well-represented by 
its OPE form. This is conventionally referred to 
as the regime of the validity of "local duality" 
(LD). As one moves to lower so, the arguments of 
Poggio, Quinn and Weinberg Q suggest that the 
OPE will break down first for those s on |s| = sq 
near the timelike real axis. We thus expect there 
to exist a regime of "intermediate" scales, sq, for 
which the OPE, while not reliable near the time- 
like real axis, continues to be reliable over most of 
the rest of the circle |s| = sq. We will refer to the 
scales for which this is true as the regime of the 
validity of "semi-local duality" (SLD). Eventu- 
ally, for sufficiently small sq, we expect the OPE 
to become unreliable over the whole of the circle 
|s| = Sq. Were all sq < rn^ to lie in this final 
regime, it would be impossible to use hadronic r 
decay data to determine parameters appearing in 
the OPE; it is therefore important to verify that 
some of the scales kinematically accessible in r 
decay lie, at the very least, in the region of valid- 
ity of SLD. 

The flavor ud V, A and V-l-A correlators pro- 
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vide an excellent laboratory for studying the na- 
ture of duality violation. The reason is that, for 
scales ~ 2 — 3 GeV^, the OPE for these correla- 
tors is completely dominated by the D = term, 
which is known to 0{al) . Since as is measured 
with good accuracy at the Z scale and the QCD 
/3-function is known to 4-loop order ||^ , so that as 
can be run down reliably to scales below 2 GeV^ , 
the OPE in these cases is known with good accu- 
racy down to scales significantly below m^. One 
can, therefore, compare the spectral and OPE in- 
tegrals for various FSER weight choices and de- 
termine the extent to which the OPE provides a 
good representation of the spectral integral data. 

The non-perturbative terms of the OPE rele- 
vant to this comparison are j7| 

11 , 



n 



(0+1) 
V/A;ud 



1 



n 



(0+1) 
V/A;ud 



D=Q 



i2g4 



[aGG) + • • • 

/ (T288 + 64) 
V 81 



(6) 



(7) 



where a = a{Q'^) = as{Q'^)/TT, p = pvsA ex- 
presses the deviation of the four-quark conden- 
sates from their vacuum saturation approxima- 
tion (VSA) values, -t- • • • in Eq. ^ stands for nu- 
merically small terms involving the quark conden- 
sates and fourth powers of the light quark masses, 
and, in Eq. (|^), the upper (lower) sign refers to 
the V (A) case. The dominant D — contri- 
bution may be expressed in terms of the Adler 
function 



(0+1) 
V/A 



D=0 



(8) 



where, in the MS scheme, Kq = Ki = 1, K2 = 
1.63982, K3 = 6.37101 [||, and the higher are 
unknown. 

To be specific, we take as OPE input (1) 
as{mr) = 0.334 ± 0.022 §,§, (2) {^G^} = 
(0.009 ± .015) GeV^[| (3) pvsA = 1 ± 5[| and use 

^This value encompasses the ranges obtamed in three 
different analyses, 0.021 ± 0.03 GeV jc|l, 0.006 ± 



0.012 GeV* Q and 0.009 ± 0.007 GeV* ||l5 
*The VSA-violating parameter, pvSA^ for the V-A corre- 



K4 = 25 ± 50 as a means of estimating the D — 
truncation error In the V-l-A case, the size of 
the VSA D — 6 contribution is quite small, so the 
conservative error on pvsA has little impact on 
the full theory errors, which are dominated by the 
uncertainty in as{mT) and the estimated trunca- 
tion error. Possible D = 8 and higher terms are 
neglected. Errors associated with the uncertainty 
in each of the inputs are combined in quadrature 
to obtain the total theory error. Spectral inte- 
grals are evaluated using the ALEPH data §,|, 
which produces the smallest error 80. The corre- 
sponding errors are computed using the ALEPH 
covariance matrix. The final errors on the ratios 
of OPE to data integrals, which provide a mea- 
sure of duality violation, are obtained by combin- 
ing theory and experimental errors in quadrature. 

Making the OPE/data comparison first for the 
V channel, we discover that sq < [0 is not 
a region of the validity of LD for the ud V corre- 
lator. Indeed, the OPE and spectral integrals for 
FESR's based on the weights w^s) = s'' (which do 
not suppress contributions from the vicinity of the 
timelike real axis) are in rather poor agreement. 
This is illustrated in Figure 1, which displays the 
OPE-to-spectral-integral-ratio (which should be 1 
if LD is valid) as a function of sq: for fc = 0, • • • , 3. 
Also shown are the corresponding results for the 
ud Y+A combination. One can see that the level 
of duality violation in the V-l-A correlator is much 
smaller than that in the V and A correlators sep- 
arately. 

The strong duality violations in the s*^- 
weighted FESR's for the ud V correlator are 
not a general feature of FESR's involving this 
correlator. Indeed, the FESR for Ry-ud, 
Eq. (|^), (whose "kinematic" weight, Wt [Ut) — 
(1 — Ur)^ (1 + 2^7-), is a linear combination of the 
four s'^-weighted FESR's shown in Figure 1) is 
known to be in good agreement with experi- 
ment [py9| . The difference has to do with the form 



lator is known to be ~ 1.6 [Q. Significant cancellation, 
however, occurs in the VSA expression for the V+A sum; 
the expanded error, ±5, on the VSA value pvSA = 1 is 
meant to deal with this in a conservative manner. 
^See also Refs. il^ for the OPAL data on the ud V and A 
spectral functions. The spectral integrals using the OPAL 
data are in good agreement with those using ALEPH data. 
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Figure 1. Ratios of OPE to data integrals as a function of sq for the flavor ud V (open circles) and V+A 
correlators (solid circles). The top left, top right, bottom left and bottom right figures correspond to 
u'(s) = 1, w{s) = s, w{s) = and w{s) = s'^, respectively. 
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of the weight: because the point s = m^, where 
the circle \s\ = crosses the timehke real axis, 
lies at the edge of hadronic phase space, Wr has 
a (double) zero at s = m^. This zero suppresses 
OPE contributions from the part of the integra- 
tion region near the timelike real axis where use 
of the OPE is expected to be potentially most 
problematic. It is presumably this feature of Wr 
which is responsible for the success of the Ry-ud 
FESR. 

If, as this argument suggests, the success of the 
Rv;ud FESR is a reflection of the localization of 
the breakdown of the OPE to the vicinity of the 
timelike real axis for scales sq — m^, then other 
FESR's with suppressions of OPE contributions 
from this region should also be well satisfied. This 
hypothesis was tested in Ref. |]l^ . FESR's for the 
V channel were considered for a range of so dif- 
ferent from m^, and a range of alternate weights, 
?«(s), still satisfying ■w{s — sq) = 0^. The result- 
ing pFESR's were all found to be well satisfied 
for all So > 2 GeV^. This suppression of duality 
violations, which is a generic feature of weights of 
the forms 

w{v) = {\-yf{l + Ay) (9) 
w{y) = {l-y){l + Ay) , (10) 

having either a single or double zero at s = sq 
(where A is a free parameter) can, in favorable 
cases (including that of the kinematic weight), 
extend to much smaller sq. This is illustrated 
for the pFESR's based on the weights w{y) = 
(1 — j/)^(l + 2y) and w{y) = 1 — y in the left 
and right panels of Fig. 2, respectively. Also 
shown for comparison in each panel are the results 
for the corresponding s*^ -weighted (unpinched) 
FESR's, where k is the degree of the pinched 
weight in question. The 'w{y) = 1 — y pFESR 
receives no contributions from the D — Q OPE 
term and hence has reduced OPE errors at low 
Sq. The absence of duality violation in this case 
extends to very low sq, despite the fact that the 
weight has only a single zero at s = sq. 

The results above provide a specific illustration 
of the dangers of attempting to extract OPE pa- 

®We refer to such weights as "pinched" weights, and the 
FESR's based on them as "pinched FESR's", or pFESR's. 
The natural variable for use in such pFESR's is j/ = s/sq. 



rameters from a sum rule analysis at intermediate 
scales when duality violations are present. Were 
one to take a typical scale sq < and "extract" 
as{so) by matching the OPE and spectral sides of 
one of the s*''-weighted FESR's, one would find a 
value which, when run up to higher scales, failed 
to agree with that obtained in direct measure- 
ments, for example, at the Z scale. In this case we 
do not actually need to know the high-scale value 
beforehand in order to expose the presence of du- 
ality violation: the SQ-dependence of the a(so) 
extracted at different sq, but using the same 
weight, can been seen to be in very poor agree- 
ment with that predicted by 4-loop running in 
QCD. Similar "internal" tests can be performed 
for any FESR. Passing such a test is of course a 
necessary, but not sufficient, condition for demon- 
strating the absence of duality violations in the 
analysis in question. 

Let us now turn to a^. In order to reduce 
the difhcult-to-quantify theoretical systematic er- 
ror associated with possible residual duality vio- 
lations, it is useful to work with correlators which 
display a reduced level of duality violation to 
begin with. The results of Figure 1 show that 
the ud V+A combination is favored in this re- 
gar dQ. Similarly, pFESR's are favored over their 
unpinched analogues at intermediate scales such 
as those forced on us by r decay kinematics. Since 
there is still significant uncertainty in the value 
of the gluon condensate, which dominates the 
non-perturbative corrections on the OPE side, 
it is useful to work with pFESR weights which 
strongly suppress the D — A contribution (Eq. ^ 
with A 2 and Eq. with A = 1). 

Bearing in mind the necessity of using the 
data not only to extract a^, but also to verify 
that duality violations are not obviously present, 
the following procedure seems optimally conser- 
vative. We first extract as using the w{y) — 
(1 - y)2(l + 2y) pFESR at sq = m^. A weight 
with a double zero and the highest possible value 
for So are chosen to increase the likelihood that 
duality violation will be negligible. Neglect of the 
£> = 8 OPE contribution (which is not suppressed 

^This combination also has reduced D = 6 contributions 
on the OPE side, and reduced errors for s > 2 GeV'^ on 
the experimental side. 
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Figiire 2. Ratios of OPE to data integrals as a function of sq for the pinched and unpinched FESR's 
involving the ud V correlator. The left figure shows the results for w{y) = (1 — y)^(l + 2y) (solid circles) 
and w{s) = (open circles), the right figure the results for w{y) = 1 — y and ■w{s) = s. 




by any factor of as for this weight choice, but 
which scales as I/sq) is also most reliable for the 
highest possible sq. Using the contour improved 
peturbation theory (CIPT) scheme [p6P] for the 
D = contribution, the result is 

a^(m^) = 0.345 ± 0.026 (11) 

where experimental and theoretical errors have 
been combined in quadrature. If instead of 
CIPT one employs fixed order perturbation the- 
ory (FOPT), the central value is reduced to 
0.326^ The method, and results, are similar to 
those of earlier analyses 0,|l|,|[|,|l|l . 

Having extracted as using only one pFESR, 
and one value of sq, we can then use other 
pFESR's and other sq values to check for pos- 
sible duality violation and/or the presence of ne- 
glected higher dimension contributionsQ. Here we 
compare the OPE predictions (obtained using the 

*FOPT involves an expansion in as{so), CIPT a summing 
of logs point-by-point around the integration contour by 
the scale choice fi^ = Q^. 

^Contributions to the OPE side from terms of dimension 
D = 2k + 2 scale as 1/sq, allowing one to distinguish 



value of as just extracted) to the corresponding 
data integrals for (1) w{y) = (1 — 2/)^(l + 2y) (at 
So < fnl) and (2) w{y) — {1 ~ y){l + y) (also as 
a function of sq)- The sq dependence of the OPE 
prediction depends crucially on the form of the 
4- loop running of a^. The results of the tests for 
possible duality violation are shown in Figure 3. 
There is no evidence of such violations, within ex- 
perimental errors, even down to rather low scales. 
While this does not prove that they are absent, 
it does gives us additional confidence in the fit 
value, Eq. (|l|). 

We conclude this section with a comment on 
the use of higher degree spectral weights for anal- 
yses of the type described above. Weights with 
degrees fc > 2 in principle involve OPE contribu- 
tions with D = 8, • • • , 2fc-l-2. Usually these contri- 
butions are assumed to be negligible, at least for 
So — "1^. Although the values of the relevant con- 
densate combinations are not known for the ud 



contamination by an operator neglected when it should 
not have been, at least if one considers the data integrals 
over a sufficiently large window of so values. 
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Figure 3. The OPE and spectral integral sides 
of two pFESR's for the ud V+A correlator. The 
value asimr) — 0.345 obtained by the fitting pro- 
cedure described in the text is used as input on 
the OPE sides. The upper curve corresponds 
to w{y) = {1 ~ y){l + y), the lower curve to 
w{y) = (l-y)2(l + 2y). 
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V+A correlator, the analogous V-A combinations 
have been extracted in Ref . . If we take these 
as representative of the scale of higher D contri- 
butions for the V+A correlator, we can make a 
rough estimate of the expected size of higher D 
contributions, and hence determine the likelihood 
that neglect of such contributions is safe. For the 
(0,0) spectral weight, w{y) = (1 — + 2y), 

this yields an estimate of ~ 0.3% of the OPE 
total for So = m^, making neglect of such contri- 
butions quite safe. For the (1,0) spectral weight, 
w{y) = (1 - y)3(l + 2y), the estimated D = 8, 10 
contributions, at so = to? i both 1.5% of the 
OPE total, while for the (2, 0) spectral weight, 
w{y) = (1 — j/)'*(l + 2y), the estimated sq — m\ 
D — 8, 10, 12 contributions are ~ 3%, 6% and 
2% of the OPE total, respectively. Because of 



the growth in the size of the relevant polynomial 
coefficients, the neglect of unknown higher D con- 
tributions thus becomes progressively less safe as 
one goes to higher spectral weights. 

3. THE STRANGE QUARK MASS 
FROM FLAVOR BREAKING IN 
HADRONIC T DECAY 

The D — (mass-independent) part of the per- 
turbative contribution to the flavor-breaking dif- 



ference An 



(J) 

V/A 



n 



(./) 

V/A;ud 



(■/) 

V/A;us 



vanishes in the SU{3) flavor limit. The for- 
mally leading (in dimension) term in the OPE 
of An^^^ is therefore the D ^ 2 term result- 
ing from flavor breaking in the m^-dependent 
perturbative contributions. To the extent that 
rh = (m„ + md)/2 « rus this contribution 
is proportional to m^. The basic idea of using 
hadronic r decay data to determine rUs is then 
to construct flavor-breaking differences of ud and 
us spectral integrals and equate these, using the 
basic FESR relation, to equivalent OPE contour 
integrals involving the unknown parameter Wg. 

A simple way to construct such flavor-breaking 
differences is to note that, from Eq. (||), the ud de- 
cay distribution, rescaled by l/|V^dP, and the us 
decay distribution, rescaled by l/|V^sp, become 
equal in the SU (3) f limit. Weighted integrals of 
the difference of the rescaled spectral functions 
then produce flavor-breaking observables of the 
type amenable to a pFESR extraction of . 

Such differences can, in principle, be con- 
structed separately for V, A and V+A combi- 
nations, and also separately for (J) = (0+1), 
(0). Since the V/A separation has not been per- 
formed for the us spectrum, and even for the 
ud spectrum the V+A spectral distribution is 
better determined than are the separate V and 
A distributions, it is preferable to work with the 
V+A combinations for both ud and us. Bearing 
in mind that aJ = 0/J = 1 spin separation is not 
currently available above ~ 1 GeV in the us chan- 
nel (where the J = spectral contributions to 
the experimental decay distributions may not be 
totally negligible) the natural choices for flavor- 
breaking observables are those constructable 
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from the kinematically- weighted combinations, 



which are determined directly from the experi- 
mental ij = ud and us decay distributions. Most 
analyses in the literature have employed the {k,n) 
spectral weight versions of this construction, 



where 



^V+A-ud ^V+A;us 



\Vud 



\Vus 



(12) 



R 



(k,n) 
V+A;ij 



7/1 \k n^^V+A;ij ,,„N 

ds{i-yr) Vr — j;'^ ■ y^^) 



The approach to determining just outlined, 
though straightforward in principle, turns out to 
have non-trivial complications. 

On the experimental side, the first complica- 
tion is the limited accuracy of the current deter- 
mination of the us spectral distribution (known 
to 6 - 8% in the K* region, and to ~ 20 - 30% 
above 1 GeV^ §). This situation will improve 
dramatically as analyses of the r decay data from 
the B factory experiments begin to come on- 
line. The second complication arises from the 
rather close cancellation between ud and us con- 
tributions occuring on the spectral sides of the 
flavor-breaking pFESR's employed in the litera- 
ture (typically to the ~ 10% or less level). Such 
close cancellation makes the ud-us spectral dif- 
ference sensitive to both uncertainties in \Vus\'^ 
and (apparently) small changes in the measured 
us branching fraction. For example, with ud-us 
cancellation to the 10% (5%) level, the 2.5% 
uncertainty in the value of |14sP @ translates 
into a 25% (50%) uncertainty in the value of the 
integrated ud-us spectral difference. The small 
differences beween the preliminary ALEPH deter- 
mination = Rv+A:us = 0.155±0.006 lH, the 
final published version, R^s = 0.161 ± 0.007 |], 
and Davier's update at Tau'2000, Rus — 0.163 ± 
0.006 111,^,1^, also produce significant shifts 
in the central value of TOs extracted in a given 
pFESR analysis for the same reason. 

Table 1 displays the central values of ob- 
tained in a number of independent analyses re- 
ported in the literature, all nominally based on 
the "same" (ALEPH) us data. Also shown 



are the central values employed for the input 
quantities Rus and \Vus\- The entries are la- 
belled by their analysis type. Table 2 shows 
the same results converted to common input, 
and also to a common D — 2 OPE truncation 
scheme, for two choices of the CKM input: (1) 
the central values of the three-family unitarity- 
constrained PDG2002 fit, |K<i| = 0.9734 and 
\Vus\ = 0.2225 (CKMU) and (2) the best fit in- 
dependent PDG2002 central fit values, \Vud\ = 
0.9749 and ^ 0.2196 (CKMN). The sensi- 

tivity to Rus and \Vus\, as well as the good con- 
sistency between the different analyses when the 
same input is employed, is evident from the ta- 
blesQ. 

On the OPE side, the major complication 
has to do with the bad behavior of the various 
weighted integrals of the D = 2 part of the longi- 
tudinal ((J) — (0)) contribution. 



dsw{y) 



An 



(0) 

V+A 



(^) 



D=2 



The expression for 



An 



(0) 

V+A 



(14) 



in terms of 



D=2 
2\ 



the running quark masses mq{Q ) and running 
couphng as(Q^), is known to ©(a^) It turns 
out that, for all of the weights studied in the liter- 
ature, the integrated version of this series is non- 
converging, even at the highest scale (sq = ml) 
allowed by kinematics [ ^2[|27| , |33[ | . This is true 
whether the integrated series is organized using 
FOPT or CIPlQ Most of the recent analyses in 
the literature deal with this problem by nonethe- 
less retaining the longitudinal contributions out 
to 0(as), combining them with the (0 + 1) con- 
tributions, and considering the sum, truncated at 
0(ag). Attempts are made to assign errors to this 
truncated sum which are sufficiently conservative 
to take into account the bad behavior of the (0) 
part of the series I2^,|9|, p8| , p2|j29[| . Exceptions are 
(1) Ref. Il^, which makes a subtraction of the 

^"^For further details on the variations in input among the 
different analyses, the impact of these variations on the 
extracted values of rris, and the conversion to common 
input and truncation scheme, see Ref. wM. 
^^As an example of the bad convergence, for sq = rri^ and 
the kinematically-weighted (0, 0) spectral weight pFESR, 
the integrated series behaves as ~ 1 -|- 0.99 + 1.24 + 1.59 
in FOPT and ~ 1 + 0.78 + 0.78 + 0.90 in GIFT [§. 
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Table 1 

Central values for TOs(2 GeV) (in MeV) in the MS scheme from various analyses. {k,0) labels analyses 
which combine the (0, 0), (1, 0) and (2, 0) spectral weights. The 1999 ALEPH analysis, labelled by (fc, n), 
is a combined analysis using the (0,0), (1,0), (2,0), (1,1), (1,2) spectral weights. The W20 label for 
KMOO denotes a non-spectral weight, the details of which may be found in Ref. the analysis in 

this case is of the non-inclusive (0 -I- 1) type. The various analyses also display some differences in their 
truncation procedures, especially in the case of ALEPH99, and some differences in their treatments of 
D = A contributions. Small differences in the central values employed for \Vud\ also exist. Details may 
be found in Ref. (2^. 



Reference Type Rus \Vus\ ms{2 GeV) [MeV] 

CKP98 13 (OTO) 0.155 0.2213 145 

ALEPH99 g {k,n) 0.161 0.2218 144 

PP99 PI (/c,0) 0.161 0.2218 114 



KKPOO 129| (0,0) 0.161 0.2218 125 

KMOO 111 W20 0.161 0.2196 115 

DHPPCOO H {k,0) 0.163 0.2225 108 



Table 2 

Central values for (2 GeV) (in MeV) obtained by updating analyses in the literature to common CKM 
input, the most recent value of Rus and, where possible, a common truncation scheme. Possible D — 8 
and higher terms have been neglected throughout. Since DHPPCOO represents a combined update of 
both ALEPH99 and PP99, we display only the update of DHPPCOO in this table. KKPOO employs an 
expansion in an effective coupling which is different from that used in the other analyses. KMOO is a 



(0 -|- 1) analysis, whereas the others are inclusive. In Ref. |30 a fc-dependent truncation scheme has been 



used for the (fc, 0) spectral weights; this scheme has been retained in converting to CKMU and CKMN 
input. A discussion of the issues underlying this choice of truncation scheme is given below. For further 
details see Ref. iGSl. 



Reference 


Type 


CKMU input 


CKMN input 


CKP98 13 


(0,0) 


116 ±31 


99 ±34 


KKPOO [H 


(0,0) 


120 ± 28 


106 ± 32 


KMOO m 


W20 


110 ± 16 


100 ± 18 


DHPPCOO H 


(0,0) 


124 ±32 


106 ± 37 




(1,0) 


113 ±32 


102 ±21 




(2,0) 


99 ±21 


91 ±21 


CDGHPPOl |§ 


(0,0) 


126 ±31 


107 ±35 




(1,0) 


116 ± 19 


106 ± 19 




(2,0) 


113 ±22 


103 ± 22 
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longitudinal contributions to the spectrum based 
on sum rule analyses of the strange scalar and 
pseudoscalar channels js^, and then works with 
pFESR's for the better behaved (0 + 1) correlator 
and (2) Ref. which truncates the (0,0) and 
(1,0) spectral weight pFESR's at 0(a?), but the 
(2,0) spectral weight pFESR at 0(a^)Q. 

The badly-converged D = 2 longitudinal OPE 
series for the {k, 0) spectral weight pFESR's turns 
out to have another problem, namely that, when 
employed with either of the truncation schemes 
noted above, it produces an unphysical decrease 
with k in the extracted value of . For such 

analyses there is thus an additional theoretical 
systematic error not accounted for in the errors 
quoted in the literature. 

It is easy to see how the problem arises. Since, 
apart from the pion pole, the ud (J) = (0) spec- 
tral function is negligible, the non-pole ("contin- 
uum") part of the spectral function of Ally^^, 

^Pv+A^ negative definite. Thus, in the (fc, 0) 
spectral weight pFESR, where it occurs weighted 
by —2yr (1 — Vt)^^^ , it yields a strictly positive 
contribution to the integrand over the whole of 
the continuum region s > sth — {mx + "t-u-)^- 
Since Q<l-yr<l - Sth/m^ = 0.87 in this 
region, the integrated continuum contribution, 
|^^('=.o)j^^ is necessarily a decreasing function of 
k, and must satisfy the rigorous inequality 



< 0.87 



(15) 



The sum of the longitudinal pion and kaon pole 
contributions is also positive, and a (slowly) de- 
creasing function of fc. 

The bound in Eq. (|l5|), though rigorous, is 
overly conservative. One would, in fact, expect 
the continuum contribution to be dominated by 
the K^{lAi{)) and i4r(1460) resonances. Since the 
masses and widths of the two resonances happen 



^■^For the (0,0) and (1,0) spectral weights, the series for 
the sum of the (0 -I- 1) and (0) D = 2 contributions is de- 
creasing with increasing order out to O(of ); for the (2, 0) 
spectral weight, however, the bad behavior of the longitu- 
dinal part wins out earlier, and the 0(af) term is larger 
than the O(o^) term. The fc-dependence of the truncation 
scheme of Ref. |BC|] for the (fc, 0) spectral weight pFESR's 
results from the ansatz of truncating the series at the point 
beyond which the terms begin to increase in size. 





~ 0.44 






L 






" ~ 0.22 

L 





to be comparable, one can make an improved es- 
timate for the relation of the [A^'^'^^j^ for differ- 
ent fc by integrating over a Breit-Wigner profile 
with the average mass and width. The result is 
that, in the limit of resonance dominance of the 
continuum longitudinal contributions, one would 
expect [ID 



(16) 



We can contrast these physical constraints with 
what is implied by the truncated OPE represen- 
tation. To do so we consider the integrated (fc, 0) 
OPE representations (using the combined fit 
central value of Ref. as input, to be spe- 
cific) and subtract from them the very accurately 
known pion and kaon pole spectral terms. This 
leaves the implicit OPE representation for the 
longitudinal continuum contributions. The re- 
sulting longitudinal contributions turn out to be 
comparable to, or larger than, the corresponding 
longitudinal pole contributions for all the cases 
under consideration (fc = 0,1,2). The longitudi- 
nal contributions are also larger than the (0 -I- 1) 
contributions for all cases. With the uniform-in- 
fc truncation scheme, one finds that the contin- 
uum longitudinal contributions predicted by the 
OPE are in the ratios 1 : 1.16 : 1.42 for the 
(0,0), (1,0) and (2,0) cases, respectively. This 
fails to satisfy even the weak rigorous inequalities 
implied by Eq. (15), let alone the expectations 
based on resonance dominance given in Eqs. (p^). 
If one instead uses the fc-dependent truncation 
scheme of Ref. ||3^, the truncated OPE implies 
continuum longitudinal contributions in the ra- 
tios 1 : 1.16 : 0.83, still rather far from those 
implied by the physical constraints. 

The failure of the OPE to satisfy the phys- 
ical constraints (which must necessarily be re- 
flected in the spectral integrals) means that the 
rris values extracted in the separate (fc, 0) spec- 
tral weight analyses will themselves have an un- 
physical fc-dependence. Since, for fixed m^, the 
higher fc OPE contributions are too large relative 
to the (0, 0) contributions, successively smaller 
values of will be required to produce a match 
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between the OPE and spectral integrals as k is 
increased. This trend is seen clearly in the re- 
sults of the uniform-in-A; truncation scheme for 
the (fc,0) analyses reported in the literature, and 
in Table 2 above. The size of the resulting sys- 
tematic uncertainty on nis (2 GeV) is likely to be 
at least 20-25 MeV for the (1, 0) and (2, 0) anal- 
yses g. 

In order to avoid the problems with the OPE 
representation of the longitudinal contributions, 
one can work with pFESR's for the (0 + 1) corre- 
lator difference alone. Since the us spin separa- 
tion above 1 GeV is not available, however, some 
external input is required in order to allow a de- 
termination of the longitudinal continuum sub- 
traction to be performed in converting from the 
measured to the purely (0-1-1) spectral distribu- 
tion. 

One possibility, tried in Ref. is to iden- 
tify experimentally the K^{U30) and 7-^(1460) 
decay modes (expected to dominate the non-pole 
part of the (J) = (0) spectral integral), and sub- 
tract these contributions, plus the longitudinal 
pole contributions, from the measured spectrum. 
Since neither the 7^0(1430) nor 7^(1460) decays 
of the T have been detected to date, this is not 
yet practical, but may well become so with the 
new B factory data. 

It is also possible to compute the continuum 
(J) = (0) subtraction if one knows the Kq{1A30) 
and isr(1460) decay constants Q]. These can be 
estimated with ^ 10 — 20% accuracy using sum 
rules for the correlators involving the divergences 
of the us vector and axial vector currents |^,^ . 
An alternate possibility for the strange scalar con- 
tinuum contributions is to use analyticity, unitar- 
ity and the existence of an Omnes relation for the 
timelike scalar Ktt form factor to compute the 
strange scalar spectral function. This was done, 
ignoring the effects of channel coupling above the 
ifg(1430), in Ref. |3^. An improved version of 
this analysis, which includes the effects of chan- 
nel coupling using a model constrained by ChPT 
and known short-distance physics, has also re- 
cently been performed The effect of chan- 
nel coupling on the spectnmr in the K^{1A30) 
region is found to be small. The 7^^0(1430) de- 
cay constants found in the sum rule analysis of 



Refs. J|34|,|35| and the coupled channel analysis of 
Ref. |37]| are in good agreement within errors, 
suggesting that the uncertainties in the theoreti- 
cal determination of the longitudinal subtraction 
are under control, certainly at the ~ 20 — 40% 
level. Since, for practical reasons, the weights 
employed in current pFESR analyses necessarily 
strongly suppress contributions from the region 
above 1 GeV^, where us spectral errors are large, 
such a level of uncertainty on the longitudinal 
subtraction produces an uncertainty in the ex- 
tracted value of rris which is negligible compared 
to the other sources of error p6| , p5|j34[ ] . 

An analysis of the (0-1-1) type was performed 
in Ref. ||2^. The analysis employs three weights, 
constructed so as to (1) strongly suppress us con- 
tributions from the region above 1 GeV^, (2) im- 
prove the convergence of the integrated (0-1-1) 
OPE D = 2 series, (3) reduce to some ex- 
tent the strong ud-us cancellation present for the 
kincmatically- weighted (0-1-1) spectral integral, 
and (4) suppress, as much as possible, poten- 
tial contributions from unknown higher dimen- 
sion [D > 6) condensates. The second point is of 
relevance because the known terms of the D = 2 
series for Ally^^^^ suggest potentially slow con- 
vergence^. Explicitly, in the MS scheme jsi]] . 



A<:i)(Q^)^^^ 

19.9332a2 - 



3m|(Q2) 
27r2Q2 



(17) 



Both the running coupling and running mass are 
known to 4- loop order |^,^ . Details of the treat- 
ment of higher D contributions may be found in 
Ref. 



1 26 



Regarding the fourth point, note that, 
for the higher degree spectral weights employed in 
the literature, the modified transverse kinematic 
weights. 



(18) 



^g;"^)(^) = y" jl- yf^' (l + 2y) , 

^''In estimating the D = 2 truncation errors it is impor- 
tant to bear in mind that, for the (A;, 0) spectral weight 
cases, there is an accidental suppression of the second or- 
der CIPT contour integral caused by cancellations between 
contributions from different parts of the contour. This 
cancellation does not persist to higher orders, so taking 
the second order term as an estimate of the truncation er- 
ror is likely to produce a significant underestimate of this 
uncertainty. 
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involve uncomfortably large coefficients of the y™ 
terms with to > 2, which correspond to integrated 
OPE contributions proportional to condensate 
combinations with D = 2m + 2 > 6. This in- 
cludes the (1, 0) and (2, 0) spectral weights, whose 
pFESR's are favorable from the point of view of 
reduced ud-us cancellation. These pFESR's are 
thus unfavorable from the point of view of pos- 
sible unknown higher dimension OPE contribu- 
tions. The weights of Ref. |^ have been designed 
to avoid large coefficients for this reason. 

Especially after the most recent update of R^^ , 
the ud-us cancellation creates sizeable errors for 
the weights wio and wiq of Ref. |2^. The re- 
maining weight, W2o{y), whose form is given ex- 
plicitly in Ref. has a profile in the spec- 
tral integral region intermediate between that of 
""^(o+i) (y) ^i^d As such, it has a lesser 
degree of ud-us cancellation than the former, but 
somewhat greater degree of cancellation than the 
latter. Whereas w[o+J) (y) = 1 — 2y — 2y'^ + 8y^ — 
7y*-\-2y^ has two uncomfortably large coefficients, 
8 and —7, however, the largest coefficient in W20, 
associated with D — 8 contributions, is 2.1. Up- 
dating this analysis for the new R^s value (fol- 
lowing the procedure described in Ref. |2|] for 
the modification of the us spectral distribution), 
truncating the D — 2 OPE series at O(a^), and 
representing the D > 6 contributions by an effec- 
tive D = 8 term, ffi to data, we find 



Figure 4. Optimized OPE ffi to the spectral inte- 
grals for the W2Q pFESR. The {D = 4) and £) = 6 
OPE contributions have been subtracted from the 
data integrals. The solid line corresponds to the 
best fit for TOj and the coefficient for the D^f / = 8 
effective operator used to represent D > 6 effects. 
The data integrals and errors correspond to the 
ALEPH data and covariance matrix. 
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TO^(2 GeV) = 123± 18± 15 MeV , (19) 
for CKMU input and 

TO, (2 GeV) = 104 ± 18 ± 17 MeV , (20) 

for CKMN input. The first error in each case is 
experimental, the second theoretical. The theo- 
retical error is dominated by the estimate of the 
error associated with truncating the D = 2 series. 
The size of the data errors, and the quality of the 
match between the OPE and hadronic sides of the 
W20 FESR which results, are shown in Figure 4. 

Possibilities for improving the (0-1-1) determi- 
nation of TOs will be described in the final section. 



4. CHIRAL LIMIT VALUES OF THE 
K ^ttt: EWP matrix ELEMENTS 

In the Standard Model, the effective 
strangeness-changing non-leptonic decay Hamil- 
tonian takes the form 

10 

= -^Vudv:, Yl + ^y^^ ' (21) 
fc=i 

where r = —VtdVil/VudV*^ and the explicit forms 
of the effective operators Qi,---,(5io may be 
found in Ref. The Wilson coefficients, Zk and 
yk, arc known to NLO in a, |||,|^,|l). The task 
of determining physical weak decay amplitudes, 
such as those relevant to e'/e, is thus reduced to 
the problem of evaluating the relevant low-energy 
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matrix elements of the operators Qk ■ 

Since CP violation in the Standard Model re- 
quires the participation of all three generations, 
e'/e receives contributions only from the glu- 
onic penguin operators Q3-6 and EWP operators 
Qj-iQ. For the large physical top quark mass, 
it turns out that the dominant contributions are 
those asociated with Qg and Qg H^- In the MS 
NDR scheme, for example ||39|, 



the 1 = 2 final tttt state, are 



e' , f ImAt 



[-2.0 GeV"^ 



x((^^)7=0|Q6|if">2 Gov(I-^^Ib) 

-0.50 GeV-3 • ((7r7r)/=2|Q8|^">2 GcV 
-0.06] , (22) 

where il/s is an isospin- breaking correction 
and At = VtdV*,. 

By good fortune it turns out to be possible 
to use hadronic r decay data to determine the 
values of the K ^ tttt EWP matrix elements 
{■kt:\Qt ^s\K) in the chiral limit. The Qs matrix el- 
ement is a crucial ingredient in our understanding 
of the Standard Model prediction for e'/e. Al- 
though the Q-j matrix is not of equivalent phe- 
nomenological interest, it is likely to be of rele- 
vance as a test of lattice techniques for the eval- 
uation of weak matrix elements, particularly the 
reliability of the extrapolations to physical light 
quark masses. We briefly review the fortuitous 
circumstances which make the evaluation of these 
matrix elements possible. 

From the explicit form of Qy, Qg, 



Qi = -^SaJ^ (1 - 75) da ^ eqqbj^ (1 - 75) qb (23) 
and 

Qs = T^Sa'j^, (1 - 75) db ^ eqqbj^ (1 - 75) qa , (24) 



where a, b are color labels, one sees that both 
transform as 8l x 8ij under SU{3)l x SU{S)r. 
As a result, the K tttt matrix elements survive 
in the chiral limit (in contrast to the gluonic pen- 
guin matrix elements) and hence can be evaluated 
using soft pion/kaon techniques. The results, for 



({nnh\Qs\K'), 
and 



■3F§ 



(Oi> + -(08) 



((7r7r)2|Q7|i^")^ = --3 (Oi)^ , 
^0 



(25) 



(26) 



where Fq is the pion decay constant in the chiral 
limit and the operators are given by Oi.g are 



0\ = '?7py'7 97^ y 9 

T-3 



-lliilb^l 97 75yg , 
Og ^ q^^X-^qq^^^r^q 

-97m75A° ^99-7^75 ^9 



(27) 



The connection to t decay data arises from the 
fact that precisely the same vacuum matrix ele- 
ments (Oi.g) determine the D = 6 part of the fla- 
vor ud V-A correlator, An„d = n[^+^^ - ^aX^J- 
Explicitly 

with 



(28) 



2 [{{2TTas^Aml 
+A,lalO,)] 



hiii) = 2 [a2(S8(Og)+Bi(Oi))] 



(29) 



where Ai^g and Bi g are coefficients depending on 
the renormalization scheme, the prescription for 
the treatment of 75, and the evanescent operator 
basis |39j . Values employing the same evanescent 
operator basis as used for computing the Wilson 
coefficients of H^,_eff ©0 and for both the HV 
and NDR prescriptions for 75, may be found in 
Ref. Q. (Oi) is suppressed relative to (Og) in 
the large Nc expansion. Dispersive determina- 
tions (to be discussed below) also indicate that 
it is much smaller numerically. Thus, both the 
dominant CP violating K tttt Qg matrix ele- 
ment and the D = 6 part of the ud V-A correlator 
are essentially determined by the same quantity, 
(Og). One can therefore construct either disper- 
sive sum rules or pFESR's for the V-A correlator 
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in an attempt to extract this matrix element us- 
ing the ud V-A hadronic r decay spectral data. 

Dispersive sum rules for (Oi^s) were first con- 
sidered in Ref. |Q and re-analyzed in Ref. ]44| , 
taking into account the D > 6 contributions dis- 
cussed in Ref. and NLO radiative corrections. 
With Ap(s) the spectral function of AHud, one 
finds, for (Oi), 

{0,),~-^{asOs),^h{fi) (30) 



where the scheme-dependent coefficient Cg is 
given in Ref. [M, and 



(31) 



with 

hifi) = ds In Ap(s) 
i7i(/i)= / dQ2Q4 [An„,(Q2)]^^^. (32) 
Similarly, for (Og), one has 

= 27ra«(M)/8(M) (33) 



where 



/8(m) 



with 



27ras(^) 



/8(A*)-ff8(A^) 



(34) 



ds — ■ — - Ap(s) 



s + Ai 

6 



D>6 



(35) 



Although the spectral integrals /i,8(m) extend 
over the whole range < s < oo and not just 
< s < m^, it turns out to be possible to use the 
relations above, together with hadronic r decay 
data, to evaluate these integrals with good accu- 
racy. This is made possible by the fact that the 
ud V-A correlator satisfies three classical chiral 
sum rules in the chiral limit, the two Weinberg 
sum rules and the sum rule for the chiral 
limit pion electromagnetic splitting These 
sum rules involve integrals over the ud V-A spec- 
tral function, Ap(s), but with weights w{s) = 1, s 



and slog(s / h?) (for any A), respectively, rather 
than those relevant to (Oi,8)- 

Any weighted integral of the form 
/q dsw{s)Ap{s) can be re-written as 



ds 



Cl + C2S - 



csslog 



dsAw{s)Ap{s) 



(36) 



where Aw{s) = w{s) — c\ — cis — c^slog{s / h?). 
The first integral in Eq. ( p6| ) is determined, for 
any ci_3, by the chiral limit values of and the 
pion squared electromagnetic mass splitting, as 
a consequence of the three classical chiral sum 
rules. So long as it is possible to make Aw{s) 



small in the region above s 



by an appro- 



priate choice of ci_3, the original integral can be 
evaluated in terms of the t spectral data and the 
chiral limit input. The details of how this "resid- 
ual weight method" (RWM) , is implemented, and 
the procedure for minimizing the resulting errors, 
may be found in Ref. El . See also Refs. Hilli 



for other approaches to employing dispersive sum 
rules, or FESR's, to evaluate the chiral limit Qt.s, 
K --^ TTTT matrix elements. 

The RWM method turns out to produce errors 
on the spectral integrals /i,8(/i) which increase as 
/J, is increased. To be able to work at lower scales, 
where the errors are small, however, one has to 
worry about the presence of the higher dimen- 
sion D > 6 contributions Hi^^[p). An alternate 
approach to determining the D — Q term in the 
ud V-A OPE is to employ appropriately designed 
pFESR's. These have the advantage of requiring 
data only over a limited range s < so, which may 
be chosen to lie entirely in the range allowed by t 
decay kinematics. In addition, they can, in princi- 
ple, be employed to extract the higher dimension 
coefficients in the V-A OPE which are needed for 
an empirical determination of -ffi,8- Such a deter- 
mination would then allow for a low-scale version 
of the RWM analysis. Since, at scales p ^ 2 GeV 
the contributions to Ii^s(fj,) associated with the 
chiral constraints are more than 50% again larger 
than the integrals of the residual weights over 
the range of s covered by the r decay data, the 
pure pFESR and low-scale RWM "hybrid" analy- 
ses are largely independent, and provide a useful 
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self-consistency check on one another. Details of 
these two analyses may be found in Refs. 
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The preliminary results from this analysis show 
excellent consistency between the two methods 
for the Qs matrix elements, and are listed in Table 
3. Note that the (Oi) contribution to ag is only 
~ —3% of that associated with (Os)- The direct 
pFESR extraction thus determines only (Os)!^- 
Thus, although one may perform a hybrid anal- 
ysis for (Oi), no pure pFESR analysis is possi- 
ble. An independent test of the hybrid extraction 
of (Oi) can, however, be obtained by performing 
a RWM analysis for (Oi) at a scale sufBciently 
high that may be taken to be zero (say 

fi ^ 4: GeV) and then evolving the results down 
to = 2 GeV using the known anomalous di- 
mension matrix . The results of this test are 
shown in Table 4. 

The results of the pFESR analysis correspond 
to an EWP contribution to e'/e, in the chiral 
limit, of 



-16.2 ±3.4) X IO- 



CS?) 



and to deviations from the vacuum saturation 
value for the Qs matrix element Bg, = 1.7 and 
1.9 for the NDR and HV schemes, respectively. 
These are a factor of ~ 2 larger than the values 
produced by most models employed previously in 
the literature. The result of Eq. ( ^7| ) is similar in 
magnitude, but opposite in sign, to the current 
experimental determination 



(18 ±4) X 10" 



(38) 



A rapid change in e'/e as TOs is varied from to 
its physical value is thus required if the Standard 
Model is to explain the experimental result. 

5. COMMENTS/PROSPECTS 

The hadronic r decay data base will increase 
significantly with the data from CLEO-C, the 
new B factory experiments and (eventually) the 
Beijing r-charm factory. In this section the ques- 
tion of what improvements in the determinations 

^*The smallness of the (Oi) contribution is a result of 
two factors: the color suppression of (Oi) relative to (Og) 
and the additional factor of as in the coefficient of (Oi). 
The results of the dispersive analysis are, in fact, used to 
evaluate the tiny (Oi) contribution to ag. 



discussed above are likely to be made possible by 
this new data is briefly discussed. 

Regarding as, it appears unlikely that signifi- 
cant further improvement can be made since the 
dominant error, even with the existing data, is 
theoretical, associated with the truncation of the 
D = series |9|jl^. This shows up in the direct 
estimates of the truncation error based on varia- 
tions in the unknown O(a^) term in the D = 
part of the Adler function and/or the renormal- 
ization scale, and also in the deviations between 
the solutions for as obtained using different meth- 
ods of handling the truncated D = series|^. 

The situation is quite different for m^, where 
significant progress can be expected in the near 
future. The B factory data should drastically re- 
duce the errors on the us spectral distribution. 
In addition to the obvious advantage of produc- 
ing a reduction in the size of the errors on rus for 
existing analyses, such improvements should also 
make possible the use of alternate sum rules not 
practically useful at present because of the size of 
the spectral integral errors when evaluated with 
current data. 

This latter point may be relevant to improv- 
ing less-than-optimal features of existing anal- 
yses. An example is the current sensitivity to 
\Vus\. At present, the large us experimental er- 
rors above the K* force one to work with weights 
which strongly suppress the high-s part of the 
spectrum. These weights turn out to produce a 
rather high level of ud-us cancellation, and hence 
a significant sensitivity to \Vus\- It is likely that, 
with improved data, the reduction in the level of 
high-s suppressions required will allow the con- 
struction of alternate weights with a reduced level 
of ud-us cancellation. 

An important point to be stressed regarding 



^^The CIPT scheme is based on the truncation of the 
Adler function at fixed order, followed by integration 
around the contour |s| = sqi the FOPT scheme on the 
expansion of the integrated Adler function to fixed order 

in as(so). Since the integral of |^Os(Q'^)j can be written 
as a series in «s(so) beginning at order k, the two schemes 
differ by terms higher order in Q!s(so). The difference be- 
tween the results Q!s(m^) = 0.345 and Os(m^) = 0.326, 
obtained from the CIPT and FOPT treatments of the 
same pFESR, at the same value of so, thus provides one 
estimate for the size of the truncation error. 
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Table 3 

pFESR and "hybrid" results for the matrix elements Afy^g = ((7r7r)/=2 1(97,8|^")2 GcV- The hybrid analysis 
is an RWM analysis at /i = 2 GeV, using pFESR input to evaluate Hi,s{2 GeV). Results are in units of 
GeV^. 



Scheme Method Ms{2 GeV) Mj{2 GeV) 



NDR 


hybrid 


1.65 ±0.45 


0.21 ±0.03 




pFESR 


1.62 ±0.34 




HV 


hybrid 


1.84 ±0.46 


0.46 ±0.08 




pFESR 


1.80 ±0.36 





Table 4 

Comparison of the "hybrid" determination of (Oi(2 GeV)) and the results obtained by evolving the 
H — 4: GeV scale RWM results to fj. — 2 GeV using the known anomalous dimension matrix. Results are 
in units of 10-4 GeV^ 



Scheme Hybrid Analysis Evolution from /i = 4 GcV 

NDR -(0.70 ±0.11) -(0.53 ±0.34) 

HV -(1.52±0.27) -(1.64±0.18) 



the extraction of rus is that, given the sizeable 
theoretical systematic problems associated with 
the behavior of the OPE representation of the lon- 
gitudinal us contributions, the inclusive spectral 
weight analysis method should almost certainly 
be abandoned. There are two ways to deal with 
the longitudinal subtraction that must be made if 
one wishes pursue sum rules based on the (0 + 1) 
correlator. 

One is the theoretical approach outlined above. 
Note that the continuum longitudinal subtrac- 
tion does not need to be known to very high ac- 
curacy since it actually has little impact on the 
(0 ± 1) sum rules, at least for the weights stud- 
ied so far in the literature. For example, using 
the K{1A6Q) and A'q(1430) decay constants ex- 
tracted from the sum rule analyses of the scalar 
and pseudoscalar channels |^ to make the sub- 
traction, one finds that the shift in the values 
of the (0+1) spectral integrals induced by the 
subtraction is less than 4% for sq in the range 
So > 2 GeV^, even for the least rapidly falling of 
the (0 + 1) weights employed in the literature (the 
(0,0) spectral weight). Thus even a 100% uncer- 



tainty in the longitudinal subtraction corresponds 
to an error vastly smaller than that associated 
with the use of the badly-behaved longitudinal 
OPE representation. 

The second approach is experimental, and 
likely to be feasible, at least to some extent, with 
the new B factory data. If current theoretical esti- 
mates for the (1430) decay constant [ ^J37| are 
correct, the r i^T-i^Q (1430) branching fraction 
should be ^ 6 X 10"^ (a factor of ~ 8 smaller than 
the current experimental upper bound), a level 
which may be reachable in the new B factory ex- 
periments. Since the Kq{1A30) decays essentially 
entirely to Kn, the only experimental complica- 
tion (apart from rate) is the 7% Ktt branch of 
the _ft'*(1410). Based on the central value of the 
ALEPH determination, B{t -* i^*(1410)i^^) = 

^1.5 ^J g^ X 10"^, one would expect a K-k back- 
ground in the K^{U30) region at the - 1 x 10"'' 
level, and with a width somewhat smaller than 
that of the K^{U30). The r VrK{U60) mode 
appears less amenable to experimental identifi- 
cation; confirmation of the theoretical estimates 
for the (1430) decay constant would, however. 
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serve as strong evidence in favor of the corre- 
sponding isr(1460) decay constant estimate. Even 
an upper bound on B{t — > Kq{1A3Q)i^t) at the 
~ 1 X 10"'' level would be useful since the impact 
of the longitudinal subtraction on the (0 + 1) nis 
extraction, at the predicted level for the two de- 
cay constants, is already known to be small. It is 
thus likely that the uncertainty in the impact of 
the longitudinal subtraction can be made rather 
small experimentally in the near future. 

Regarding the K rnr matrix elements of 
the EWP operators, the main improvement to be 
hoped for is reduced errors in the ud V-A spectral 
function in the region above 2 GeV^ This would 
require an improved V/A separation of the con- 
tributions to states containing a KK pair. Such 
an improvement might be possible with the en- 
hanced statistics expected from the B factory ex- 
periments, at least for quasi-two-body modes. A 
reduction of the errors on the V-A spectral func- 
tion would allow one not only to reduce the errors 
on the extracted EWP matrix elements, but also 
to sharpen tests for the absence of duality viola- 
tion. Recall that both the ud V and A correlators 
display, in general, significant duality violation for 
So < m^. Since the V and A duality violations 
largely cancel in the V-l-A sum, the intrinsic level 
of duality violation in the V-A correlator must 
be expected to be large. The analysis of Ref. 
tests for the presence of such violations in the 
pFESR's employed, and finds none, within ex- 
perimental errors, for So above ~ 1.8 GeV^ The 
larger errors in the upper part of the so window, 
however, leave room for improvement in these 
tests, particularly in the case of the pFESR's em- 
ployed to extract the D > 6 contributions needed 
as input to the hybrid RWM/pFESR version of 
the analysis. 

I would like conclude by acknowledging use- 
ful exchanges with Andreas Hocker and Shaomin 
Chen concerning the ALEPH data tabulation, 
and especially the work of my collaborators on 
some of the topics reported here (Joachim Kam- 
bor on the light quark masses and related is- 
sues, and Vincenzo Cirigliano, John Donoghue 
and Gene Golowich of the K ^ tttt matrix el- 
ements of the EWP operators). 
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